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5 FIELD OF THE INVENTION 

The present invention is directed to a method of adaptive digital 
filtering and more particularly to such a method utilizing a reduced rank 
technique. 

10 BACKGROUND OF THE INVENTION 

Adaptive digital filters are known for use in various signal 
processing applications including speech and radar processing, adaptive 
beamforming, echo cancellation and equalization. These filters estimate 

15 the optimal filter coefficients from observed data, i.e. data representing 

the received signal to be processed. 

Reduced rank filters estimate the filter coefficients with a 
relatively small amount of observed data. Various reduced rank 
techniques are known including Principal Components, Cross-Spectral 

20 and Partial Despreading methods. The former two methods require an 

explicit estimate of the si^al subspace via an eigen-decomposition of 
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the input covariance matrix which is extremely complex. Although the 
latter technique is much simpler, it does not achieve near full rank 
performance when the filter rank D is significantly less than the full rank 
N. 

Other known reduced rank techniques include a multi-stage 
Wiener filter as proposed by Goldstein et al. in "A Multistage 
Representation Of The Wiener Filter Based On Orthogonal Projections" 
IEEE Trans. Inform. Theory, 44 (7), November, 1998 and an adaptive 
interference suppression algorithm as proposed by Honig et al. in 
"Adaptive Reduced-Rank Residual Correlation Algorithms For DS- 
CDMA Interference Suppression" In Proc. 32 Asilomar Conf Signals, 
Syst. Comput., November, 1998. These methods perform well and do 
not require an eigen-decomposition. The present invention is an 
improvement of these latter two methods. 

BRIEF SUMMARY OF THE INVENTION 

In accordance with the present invention, the disadvantages of 
prior adaptive digital filters noted above have been overcome. The 
reduced rank adaptive filter of the present invention projects a received 
signal onto a lower dimensional subspace where the reduced rank 
subspace is iteratively constructed by multiplying the last basis vector by 
the received sample covariance matrix. The reduced rank filter of the 
present invention achieves full rank performance with a projected signal 
subspace having a significantly smaller dimension than the dimension of 
the received signal subspace. Moreover, for interference suppression 
applications, the optimum filter rank does not substantially increase with 
the dimension of the received signal subspace. As such, the method of 
the present invention enables faster tracking and convergence with 
significantly less training samples than can be achieved with prior 
techniques. 
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More particularly, the method of the present invention filters 
successive, received signal samples to provide a desired signal wherein a 
group of N successive samples form a received sample vector of digital 
data having a N x 1 dimension. The method includes generating D + 1 
basis vectors where D is less than N and represents the reduced rank 
dimension. Each successive basis vector is generated from a given or an 
estimated steering vector and successively greater powers of a 
CO variance matrix for a sequence of received sample vectors of data. 
The first basis vector is formed from the steering vector. A reduced rank 
vector of digital data having a D x 1 dimension is generated from a 
matrix of D basis vectors and the received sample vector of data. D 
filter coefficients are generated from correlations between pairs of the 
basis vectors. The desired signal output from the filter is generated from 
the filter coefficients and the reduced rank vector of data. 

The method of the present invention can be used in any adaptive 
filtering application including echo and noise cancellation, channel 
equalization, radar processing, adaptive antenna arrays and interference 
suppression. The method of the present invention is particularly 
advantageous for applications where long filter lengths are required and 
fast convergence and tracking are important. One such application is for 
space-time interference suppression in wireless Code-Division Multiple 
Access (CDMA) systems. 

These and other advantages and novel features of the present 
invention, as well as details of an illustrated embodiment thereof, will be 
more fully understood from the following description and drawings. 

BRffiF DESCRIPTION OF THE SEVERAL VIEWS 
OF THE DRAWING 

Fig. 1 is a block diagram of the reduced rank adaptive filter of 
the present invention for a wireless CDMA communication system; 
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Fig. 2 is a flow chart illustrating the steps of the reduced rank 
adaptive digital filtering technique of Fig. 1; 

Fig. 3 is a block diagram illustrating the generation of the 
received signal sample covariance matrix and steering vector; 
5 Fig. 4 is a diagram illustrating the generation of the basis vectors; 

and 

Fig. 5 is a flow chart illustrating the method of selecting the 
dimension of the reduced rank filter. 

1 0 DETAILED DESCRIPTION OF THE INVENTION 

The reduced rank adaptive digital filtering technique 10 of the 
present invention is illustrated in Fig. 1 for interference suppression in a 
wireless CDMA communication system having a receiver 12. Although 

15 illustrated for a wireless communication system, the filtering technique 

of the present invention can be used in any adaptive digital filter 
application. The receiver 12 picks up a received signal via an antenna 
14 coupled to a radio frequency (RF) amplifier and mixer 16. The 
output of the amplifier/mixer 16 is coupled to an analog to digital 

20 converter 18 to provide a serial stream of baseband digital data samples 

representing the received signal. The output of the digital to analog 
converter 18 is coupled to a serial to parallel converter 20. The serial to 
parallel converter 20 takes N serial samples and outputs the N samples 
in parallel to provide a received sample vector r(0 of digital data 

25 having a dimension of N x 1 . 

The received sample vector r(/), which is referred to as an 
observation signal, is coupled to an adaptive digital filter 22, the output 
of which is an estimate of the desired signal 6(0 > referred to as the 

approximate desired signal bii) . A generator 24 produces a matrix M of 
30 basis vectors with a N x D dimension where D is less than N and 

represents the reduced rank dimension. The generator 24 also produces 
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a vector of filter coefficients c of dimension D x 1 from a sequence 
received sample vectors r(i) for / = 1, i? and from a sequence of the 
decision derived desired signals b(i), i = 1, or a sequence of given 
desired signals or training symbols Z?(/), / = 1, .. 5 . The given training 
5 symbols b(i) initially used by the generator 24 may be stored, for 

example, in a receiver memory (not shown). 

The adaptive filter 22 generates a reduced rank or projected 
received signal vector r (/) of dimension D x 1 by multiplying the 
Hermitian transpose of the matrix M by the received sample vector r(i) 
10 as follows. 

r(z) = MV(0 

The adaptive filter 22 generates the approximate desired signal b(i) by 
multiplying the Hermitian transpose of the filter coefficient vector by the 
reduced rank vector as follows. 

15 ki) = c'r(i) 

The approximate desired signal b(i) is applied to a conventional sheer 
23 that essentially rounds the approximate desired signal to a 

desired signal level b(i) which may be used to generate the filter 
coefficients in a decision directed training mode. 

20 The generator 24 generates D + 1 basis vectors where each 

successive basis vector is generated from a given or an estimated 
steering vector p and successively greater powers of a sample covariance 
matrix R for a sequence of B received sample vectors r(/), / = 1, B . 
The first basis vector is formed from the steering vector p as described 

25 in detail below with reference to Fig. 4. The generator 24 generates a 

D X 1 vector of filter coefficients c from a matrix and vector of 
correlations between pairs of generated basis vectors. 
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The generator 24 and adaptive digital filter 22 operate in 
accordance with the method depicted in Fig. 2. At block 26, the 
generator 24 determines the sample covariance matrix R for the 
sequence of received sample vectors r(£), ^ = 1, ^ and the steering 
5 vector p from the sequence of received sample vectors r(0, / = 1, . . . , 5> 

and the sequence of training/desired signals / = 1, ... in 

accordance v^ith the block diagram of Fig. 3. Specifically, at block 28, 
the Hermitian transpose of the received sample vector is formed and at 
30 is multiplied by the received sample vector and a weighting factor of 

10 1 - w. The resulting product is added at 32 to the product formed at 34 

by multiplying the weighting factor w by the previous covariance matrix 
so as to form the sample covariance matrix /?- for the sequence of 
received sample vectors r(l), r(i) . At 36, the conjugate of the given 
or desired signal b(i) is formed and at 38 is multiplied by the received 

15 sample vector r(i) and the weighting factor I - w. The resulting 

product is added at 40 to the product of the previous steering vector p.^^ 
and the weighting factor w to form the new steering vector /?. . The 
steering vector generated from the sequence of given training 
symbols / - 1, S or from the sequence of desired filter outputs 

20 = 1, 5 as shown in Fig. 3 is an example of an "estimated" 

steering vector. It should be appreciated that Fig. 3 is only one example 
of a technique for generating the estimated steering vector. Other 
known techniques for generating an estimated steering vector can be 
used as well. A "given" steering vector can be used in accordance with 

25 the present invention as well. A "given" steering vector is a steering 

vector that is already known to the receiver and may be stored, for 
example, in a memory or buffer of the receiver. 
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At block 44, the generator 24 generates D + 1 basis vectors Vq 
through Vjr,. The basis vectors are generated in accordance with the 
diagram depicted in Fig. 4. The first or initial basis vector is set 
equal to the steering vector p. Each successive basis vector is formed by 
multiplying the immediately preceding basis vector by the covariance 
matrix R so that the basis vectors are generated as follows. 

The generator 24 at block 46 forms the matrix M of basis vectors from 
the basis vectors Vq through v^,! so that the matrix M has a NxD 
dimension. Thereafter, the generator 24 determines the correlation 
scalars through y^o-x from the basis vectors Vq through 
determined at block 44. The correlation scalar y^^j is generated by 
multiplying the Hermitian transpose of the basis vector v. by the basis 

vector Vj as follows. 

^ t t t t t 

Also = Vq and Vq is assumed to be normalized to that = 1 • 

The generator 24 at block 54 forms a correlation matrix T from 
the correlation scalars y^ through /jo-i follows. 

p _ ^2 Yz *''Yd+\ 

_Yd Yd+I' "Y2D-i_ 

At block 54 the generator 24 also generates a reduced rank correlation 
vector y which is a D x 1 vector formed of the correlation scalars 

through y^_i as follows. 
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The generator 24 at block 56 then generates the filter coefficient vector 
c which is a D X 1 vector generated by solving the following set of 
linear equations. 
5 Tc ^ y 

The matrix M of basis vectors formed at block 46 and the filter 
coefficients c generated at block 56 are applied to the adaptive digital 
filter 22 so as to generate the reduced rank vector r{i) and the 

approximate desired signal b{i) output from the filter 22 in accordance 
10 with blocks 58 and 60. In particular, at block 58, the adaptive digital 

filter 22 generates the reduced rank vector r(/) by multiplying the 
Hermitian transpose of the matrix M of basis vectors by the received 
sample vector of digital data r(z) as follows. 

r(0 = MV(0 

15 At block 60, the filter 22 generates the approximate desired signal b{i) 

by multiplying the Hermitian transpose of the D x 1 filter coefficient 
vector c by the reduced rank vector r{i) as follows. 

b{i) = c'r(i) 

It has been found that the reduced rank adaptive digital filter 
20 described above achieves near full rank performance with a value of D 

equal to 8 or less for a Direct Sequence (DS) Code Division Multiple 
Access (CDMA) communication system. For other applications, the 
value of D can be determined in accordance with the method depicted in 
Fig. 5. As shown in Fig. 5, at block 60 an index n is initialized to 1. 
25 Thereafter, at block 62 the /i"" basis vector v„ is generated by 
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multiplying the previous basis vector v„_i by the covariance matrix R 
where the first basis vector is the steering vector p so that 

v„ = R^p 

as discussed above. Thereafter, a matrix M^_^ of basis vectors is formed 
at block 64 as follows. 

At block 66, the orthogonal projection m„ , which is a N x 1 vector, is 
generated as follows. 

"„ = v„ -M^_,{M\-xM,_,rMjv^ 
The respective lengths 1 and g of the vectors w„ and v„ are computed at 
block 67 in accordance with the following equations. 

' = kll 

s = IK II 

where 1 and g are scalars. Thereafter, at block 68, it is determined 
whether the length 1 divided by the length g is greater than a small 
constant S , where S may be set equal to 0.01 for example. If so, n is 
incremented by 1 and steps 62 through 68 are repeated. Steps 62 
through 70 are repeated until 1 divided by g is not greater than S for a 
given value of n. At that point, D is set equal to n - 1. 

In an alternative method, the value of D is chosen to minimize 
the "a posteriori Least Squares cost function C" where 

C = Y.Hm)^c\r{mf 

where is the reduced rank filter coefficient vector generated from the 
sequence of b(l), . . . b{B) and the sequence of r(l), . . . r(B) and < . 

The reduced rank adaptive filter of the present invention projects 
the N X 1 received signal vector r(i) onto the lower dimensional 
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subspace which results in a D x 1 vector r (z) where D can be much less 
than N while still allowing the filter 10 to achieve near full rank 
performance. For interference suppression applications, the optimum 
filter rank D does not substantially increase with the dimension of the 
5 signal subspace N. This enables much faster tracking and convergence 

than can be achieved with prior techniques. Furthermore, the method of 
the present invention does not rely on an estimate of the signal space via, 
for example, eigen-decomposition so that the complexity of the present 
method is substantially reduced. Although the method was described in 

10 detail for a wireless communication system application, it can be used in 

any adaptive filtering application including, echo and noise cancellation, 
channel equalization, radar processing, adaptive antennae arrays, 
interference suppression, etc. 

Many modifications and variations of the present invention are 

15 possible in light of the above teachings. Thus, it is to be understood 

that, within the scope of the appended claims, the invention may be 
practiced otherwise than as described hereinabove. 

What is claimed and desired to be secured by Letters Patent is: 



20 



